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ABSTRACT 

We study the implication of the collapsar model for Long Gamma-Ray Bursts (LGRBs) 
on the metallicity properties of the host galaxies, by combining high-resolution N-body sim- 
ulations with semi-analytic models of galaxy formation. The cosmological model that we use 
reproduces the Fundamental Metallicity Relation recently discovered for the SDSS galaxies, 
whereby the metallicity decreases with increasing Star Formation Rate for galaxies of a given 
stellar mass. We select host galaxies housing pockets of gas-particles, young and with different 
thresholds in metallicities, that can be sites of LRGB events, according to the collapsar model. 
The simulated samples are compared with 18 observed LGRB hosts in the aim at discriminat- 
ing whether the metallicity is a primary parameter We find that a threshold in metallicity for 
the LGRB progenitors, within the model galaxies, is not necessary in order to reproduce the 
observed distribution of host metallicities. The low metallicities of observed LGRB hosts is 
a consequence of the high star formation environment. The star formation rate appears to be 
the primary parameter to generate a burst event. Finally, we show that only a few LGRBs are 
observed in massive, highly extincted galaxies, while these galaxies are expected to produce 
many such events. We identify these missing events with the fraction of dark LGRBs. 

Key words: gamma-rays: bursts - host galaxies . 



1 INTRODUCTION 

The study of the host galaxies of Gamma-Ray Bursts (GRBs) is 
essential for understanding their nature. Current observations re- 
veal that long-duration GRBs occur in star-forming galaxies, in 
contrast to short-duration GRBs that are found in both early- and 
late-type galaxies. The star formation rate (SFR) in th e short GRB 
host g alaxies is often lower than in the LGRB's hosts jBerger et al.l 
l200d [ and references therein). The difference in the observed host 
properties for short and long GRBs supports the idea that short 
and long GRBs have different progenitors. Long-duration GRBs 
are bel ieved to arise from the death of massive stars (collapsar 
model: [ Conselice et al ]|2005l: iFruchter et al.ll2006t iTanvir & Levaril 
l2007l : IWainwright et al.ll2'o07l and references therein), most likely 
Wolf-Rayet stars, while short-duration GRBs are likely produced 
by th e coalescence of compa ct objects in binaries jLi & Paczvnskil 
Il998l : lO'Shaughnessv et aUlioOS) . The discovery of the connec- 
tion between long G R Bs and core-collapse Type Ibc supemovae 
jGalama et al.ll 19981 : 10 120061 : IWooslev & Hege3l2006l and refer- 
ences therein) supports the collapsar model of long-duration GRBs 
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jMacFadven & Wooslevl[T99^ : iMacFadven et alj|200lh . So far no 
supemovae have been found to be associated with short-duration 
GRBs. 

In current core-collapse models for long-duration GRBs 
(LGRBs hereafter), young stars with initial mass > 30 Mq 
should be able to create a black hole (BH) remnant. If the 
collapsing core has high angular momentum, the formation of 
the BH may be accompanied by a LGRB event (Woosley 
1993; MacFadyen & Woosley 1999). The wind-driven mass 
loss in massive stars is a function of the metal content: high 
metallicity stars have strong stellar winds with main losses 
of both mass and angular mom entum. Instead, when m etal- 
licities are low (0.1 - 0.3 Zq; IWooslev & HegeJ l200d ). the 
specific angular momentum of the progenitor star allows for 
the loss of the hydrogen envelope w hile preserving the he- 
lium core that can still carry rotation ( Woosley & HegeJ 12006 : 

' anl2 



Fryer. Woosley & Hartmamilll999l:lYoon. Langer & Normanl2006 



Yoon. Langer. Ca ntiello. Woosley & Glatzmaieil2008 ). The loss of 
the hydrogen envelope reduces the material that the jet needs to 
cross in order to escape, while the heliu m core should be massive 
enough to collapse and power a LGRB. IVink & de KoteJ j2005h 
explored LGRB models with even lower metallicites {Z/Zq ^ 
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10^^), showing that they have too low mass-loss rates to power 
a LGRB. Thus, it may be possible to produce a LGRB in the metal- 
licity range around (0.1 — 0.3) of the solar metallicity. 

The relationship between LGRB progenitors and their host 
environment has become matter of a hot debate, in recent years. 
In particular an open question concerns the metallicity of the 
LGRB host galaxies: is the metallicity a discriminant for the 
formation of a LGRB? Do LGRBs form preferentially in metal 
poorer galaxies? The observational information gathered so far in- 
dicates that most LGRBs are found in faint, star forming galax- 
ies dominated by young stellar populations with sub-solar gas- 
phase metallicities, although there are a few host galaxies with 
highe r metal conten t (Prochaska et al. 2004; Wolf & Podsiadl owsk i 



20071: lFYnbo etalJl2 006: Pri ce et al.l l200l ISavaglio et al.i l2003l : 
SavaglidboodlStanek et alji2006l : lLevesaue et al.ll20ld and refer- 
ences therein). We have to stress that low-metallicity progenitors do 
not necessary imply low-metallicity host galaxies. Inde ed, owing to 
the ex istence of metallicity gradients inside galaxies jArtale et al.l 
l201lh . LGRBs could form from low-metallicity progenitor stars 
collapsing on metal de ficient clouds also in hosts with relatively 
high mean metallicities dCampisi et a"Lll2009l) . 

The existence of a possible metallicity discriminant in the 
GRB selected galaxies with respect to the overall field galaxy 
population can be tested comparing observed phenomenological 
relations in the two samples. In particular, the relation between 
stellar mass M* and metallicity Z provides a good des cription 
of the properties of nearby galaxies dTremonti et al.ll2004l) . It has 
been shown that s i milar relations hold a lso at higher redshifts 
I Jabran Zahid et al.l 1201 d: lErb et al.l |2006|) up to 



I Mannucci et alj|200^ Maiolino et all 200 Sh with an evolution to 



wards lower metallicities with increasing redshift. 

Recently, iHan et all ( I2OI0I ) and iLevesgue et alj feOlOl ) have 
compared the M* — Z relation for LGRB host galaxies with sam- 
ples from the Sloan Digital Sky Survey (S PSS) representative o f 
the general star-forming galaxy population ('Tremonti et al."2004|). 
They find that the metallicity content of low redshift LGRB hosts 
tends to fall off the A/, — Z relation, and suggest that LGRBs 
do occur preferentially in lower metallici ty galaxies. In order to 
further explore the origin of this offset, iMannucci et al.l ( 1201 ll) 
compared the observed properties of LGRB hosts with those of 
field g alaxies in light of the n ew Fundamental Metallicity Relation 
(FMR: lMannucci et alj2oTol) . FMR is a tight relation between stel- 
lar mass Af«, metallicity Z and SFR. Local SDSS galaxies define 
a surface in the 3D space of these three quantities, with metallic- 
ity well determined by the stellar mass a nd SFR (with a scatte r 
around this surface of about 0.05 dex; see IMannucci et"aLll2010l) . 
The key parameter that tighten the relation is a linear combina- 
tion (in log scale) of stellar mass and star formation rate: /10.32 = 
logAf, - 0.32 log(SFR). Mannucci et al. (2011) compared the 
FMR for galaxies with stellar masses down to 10*'^ Mq with the 
metallicity properties of 18 host galaxies of LGRBs with redshift 
2 < 1, for which Af*, SFR and Z were known. They found that 
while LGRB host galaxies show a systematic offset toward lower 
metallicities with respect to the A/, — Z relation of field galaxies, 
no offset is present on the FMR. This indicates that the deviation 
relative to Af« — Z is due to the higher-than-average SFR observed 
in the hosts of the GRBs: the lower metallicity content is a con- 
sequence of the occurrence of LGRBs in low-mass, actively star- 
foiTning galaxies which are, per se, metal poorer. In other words, 
the metallicity observed in LGRB hosts is exactly what is expected 
on the basis of their mass and SFR, with no apparen t bias toward 
lower metallicities. In addition. lMannucci et al.l ( l201 ih confirm that 



LGRB hosts have in general higher specific SFRs, i.e. higher SFR 
per unit stellar mass (SSFR hereon) as in ISavagUo et al] (|2009|) . 
suggesting that the condition for a galaxy to host a LGRB could be 
related to its ability to form stars in a efficient way. 

Focus of this paper is at investigating whether the low metal- 
licities o£jhe_2rogenitorsJoLGRBs requested by the collapsar 
model dWooslev & Hegej2006l) are in agreement with the observed 
properties of the galaxy hosts, in light of the findings by Man- 
nucci et al. (2010). The questions to answer are: do LGRBs pref- 
erentially select hosts with mean metallicities lower-than-average 
among star forming galaxies? Is the SSFR the primary physical 
parameter for the formation of a LGRB rather than the metal con- 
tent? To this aim we use a simulated catalogue of galaxies, con- 
structed by combining high-resolution N-body simulations with a 
semi-analytic prescription of galaxy formation. This allows us to 
identify different galaxies and to select those housing pockets of 
gas clouds with a metallicity threshold as requested by the collap- 
sar model Wooslev & Heger (2006). 

The paper is organized as follows: in Section 2, we test the 
simulation outputs against the observed FMR and Af* — Z relation, 
and show that the simulated galaxies provide a good description of 
the two relations. Then in Section |3] we analyze the properties of 
LGRB host galaxies, defined in the context of the collapsar model, 
and suitably selected within the simulated cosmic volume. In Sec- 
tion 4, we explore the loci in the SSFR— Z and /io.32 — Z planes 
where we have the higher probability to host a LGRB. Finally, in 
Section|5]we present our conclusions. 



2 THE FMR OF SIMULATED GALAXY CATALOGUES 

In this Section we test whether the simulated galaxies reproduce 
the FMR. Read e rs inte r ested to the ga l axy's c a talogues are referred 
to IWang et al.l ( |2008|) . ICroton et al.l ( l2006h . iDe Lucia & BlaizotI 
and references therein for details on the physical processes 
explicitly modelled. Here, we give a short description of the cata- 
logue, focusing on the treatment of metal enrichment. 

The galaxy catalogue used was constructed by IWang et 

( I2OO8I) with cos mological parameter s consistent with the third-year 
WMAP results dSpergel et af]|2007h . The simulation corresponds 
to a box of 125fe~^Mpc comoving length and a particle mass 
7.78 X 10** Mq. Simulation data were stored in 64 outputs and 
analysed with the post-processing software originally developed for 
the Millennium Simulation (Springel et al. 2005). Merging history 
trees for self-bound structures extracted from the simulations were 
used as input for the Munich semi-ana l ytic m odel of galaxy forma- 
tion described in iDe Lucia & BlaizotI aioOTh . In particular, in the 
simulated galaxies the metals are produced primarily in the super- 
novae and are deposited directly in the cold gas prese nt in the disc 
of the galaxy (istantaneous recycling approximation; ICroton et al.l 
I2OO6I; iDe Lucia etani2004l) .FI 

We limit our analysis to galaxies with stellar mass larger than 
8 X 10* M0, which is above the resolution limit of the N-body sim- 
ulations used. In order to be consistent with the observed local sam- 
ple of SDSS galaxies, we select, in our simulated c osmological vol- 
ume, galaxies with redshift 0.07 < z ^ 0.3 as in IMannucci et al.l 
(2010), for which the SFR, stellar mass and mean metallicity are 
known. The resulting sample includes 1075878 galaxies. 



^ The interpretation of the FMR in the context of the simulated catalogues 
is beyond the goal of present work. 
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Figure 1. Simulated galaxies in the local Universe (0.07 < z 0.3) (lines) and galaxies of the SDSS (points in panels a,b,c); lines/points refer to the median 
values. Panel (a) shows the metallicity Z versus stellar mass M, (in units of M©); colored lines/points show the Af* — Z relation, in three selected SFR bins. 
Solid line is the observed relation (eq. 1 of the footnote). Panel (b) shows the metallicity versus SFR (in units of Mq yr~^), for galaxies in three different 
Alt, bins; Panel (c) shows the Fundamental Metallicity Relation, i.e. metallicity Z versus /^o.32 = log(M*) — 0.32 log(SFR) in three different SFR bins. 
Dark-solid line refers to the FMR relation from Mannucci et al. (201 1 ; eq. 2 of the footnote). 



Fig[T] shows the result of our comparison: panel (a) refers to 
the M, — Z plane where the metal content Z is express in terms 
of the oxygen to hydrogen abundance ratio Z — 12 + log(O/H)0 
(left 2/— axis), and in solar units (right y— axis); panel (b) refers to 
the SFR— ^ plane; and in panel (c) refers to the /io.32 — Z plane 
where the FMR is defined □; In the first three panels, lines refer to 

^ The best fit to the observed M , — Z relation for galaxies in the local 
universe is jMannucci et al.l2O10l) : 

12 + log(0/H) = 8. 96 + 0. 31m-0.23m2-0.017m^ + 0. 04601" (1) 

where m = [log( A/, /Mq ) — 10] . We will refer to eq. 1 hereon as Af* — Z 
relation. 

^ The analytical expression of the FMR for local galaxies is 



the simulated galaxies and points to SDSS galaxies. In more detail, 
panel (a) shows the median metallicity for the simulated galaxies in 
different SFR bins as a function of the stellar mass. On average each 
bin contains more than 10'^ galaxies. Galaxies with higher SFRs 
show systematically lower metallicities at a fixed A/*. Indeed this 

iMannucci et alj201lh : 

12 + log(0/H) = 8.90 + 0.37m - 0.14s - 0.19 

+0.12ms - 0.054 for /^o.32 > 9.5 

= 8.93 + 0.51(^io.32 - 10) for /xo.32 < 9.5 

(2) 

where m = [log(Af*/M0)- 10] , and s = log(SFR) with SFR expressed 
in units of Mq yr~ ^ . 
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trend is clearer in panel (b) where the SFR— Z relation is plotted 
in three different Af* bins. This finding is in agreement with the 
behavior of the observed SDSS galaxies in the same redshift range 
(depicted as points). As shown in panel (c), once the dependency on 
the SFR and A/, is included into the parameter the medians, 

in each SFR bin, collapse into a single, well defined line consistent 
with the observed FMR shown as dark-solid line. 

In conclusion our tests have shown that the simulated galaxy 
population provides a good description of the observed properties 
of local SDSS field galaxies. 



3 SIMULATED LGRB HOST GALAXIES 

In this Section, we use the simulated galaxy sample to select possi- 
ble LGRB host s in the light of the collap sar scenario. 

Following IWooslev & Hege3 l l2006l) . LGRBs occur more fa- 
vorably in stellar environments with metallicities lower than 0.1- 
0.3 Zq. Accordingly, we extracted from the simulated catalogue 
the galaxies which contain pockets of gas "particles" that can be 
sites to progenitors to LGRB events. Pockets have a metallicity 
Zprog and age tagc, and the host has metallicity ^host computed 
as mean of the metallicity of all pockets presen t inside the simu- 
lated galaxy. According to ICampisi et all ( l2009h , we create three 
host's catalogues with different properties: 

(i) HOSTl, obtained by selecting galaxies containing stars with age 
tage < t, = 5 X 10'' yr; 

(ii) HOST2, including galaxies with stars of age tagc < t* and 
metallicity Zprog ^ 0.3^0; 

(iii) HOST3, defined by selecting galaxies containing stars with age 
tage < t* and metallicity Zprog ^ 0-1 Zq. 

We computed the number of stars ending their lives as LGRBs, 
assuming a Salpetejf] Initial Mass Ftmction (IMF), and on aver- 
age (o ver all cosmic times) one LGRB event every 1000 super- 
novae jPorciani & MadajlioOll : ICampisi et al.ll20Ilh . The num- 
ber of expected LGRBs in the fc-th galaxy can be computed as 
Nk,i <x ("bh,! Mf^i^k{z), where (BH,i is the fraction of stellar 
mass that will produce black holes in case i (where i runs over the 
three catalogues, i = 1, 2, 3) and Af,.i_fc is the stellar mass in the 
fc-th galaxy t hat satisfies our se lection criteria for the i-subsample 
(see Sec. 3 in ICampisi et ai]|2009i for details). We point out that we 
select host galaxy candidates assuming constraints on the metallic- 
ity of LGRB progenitor stars, not on the overall metallicity Zhost 
of the host galaxy. 

3.1 M^- Z and FMR of simulated LGRB hosts 

In this Section, we focus on the position of LGRB hosts in the 
planes shown in Fig.l and on the comparison with the observed 
LGRB hosts. 

We first include the simulated galaxies belonging to the three 
samples HOSTl-2-3 in the Af* — Z and fio.32 — Z planes. The 
metallicity, in this case, refers to the galaxy as a whole, and is ex- 
tracted as Zhost from the catalogue. For each sample, we computed 
the median metallicity values as a function of the stellar mass A/, 
and similarly, for FMR, as a function of /io.32. In each sample, the 
median is obtained by weighting the host galaxies with the rate of 
occurrence of LGRBs in their pockets. In addition we restrict to 

* We check that an alternative IMF, as the Chabrier, does not change the 
selection of host galaxies in our simulation. 



galaxies al z ~ 0.5 being the mean redshift of the observed sample 
of LGRB hosts. 

Fig. [2] shows the results for HOSTl-2-3 and the comparison 
with the LGRB's data points. The left (right) panel refers to the 
Af« — Z relation (FMR), and lines refer to the median computed 
for each sample as described in the figure caption. The black line 
in the right (left) panel refers to the Mt, — Z re lation observed 
for galaxies at z ~ 0.8 Ijabran Zahid et al1l2010l (to FMR as in 
iMaimucci et al.ll201ll) . The blue data points and blue dotted-line 
refer to the observed sample of 18 LGRB hosts. 

In the left panel of Fig. 2, HOSTl-2-3 show a systematic offset 
toward lower metallicities compared to the observed Ad, — Z rela- 
tion for field galaxies. H0ST2 and HOST3, corresponding to low 
cut-off in Zprog, display a flat A/, — Z relation. In these two sam- 
ples, single simulated LGRB host galaxies with large M, may have 
high metallicities (i.e. high Zi^ost ) as illustrated and discussed in the 
Appendix A. However, their contribution to the median, weighted 
for the rate of LGRB housed in the galaxy, is negligible as they are 
found to be site of only few LGRB events. In the right panel of Fig. 
2, HOSTl-2-3 are compared with the FMR. HOSTl reproduces the 
observed relation for field galaxies quite well. By contrast H0ST2 
and H0ST3 show significant deviations from the FMR, as in the 
Af, - Z plane. 

The data of the 18 LGRB hosts are plotted in Fig. 2, and the 
blue-dotted line describes the best linear fit to the data (blue-poin ts 
with associated errors are as reported in iMannucci et al.l (I2OIII) ). 
In the M, - Z plane, GRB060218 is included in the fit though it 
lies at log(A/«/MQ) — 7.78 which is outside the range of mass 
accessible from the simulated volume. 

In both planes, HOSTl is the sample providing the best de- 
scription of the available data. In particular, this sample reproduces 
well the FMR for the observed LGRB hosts and the systematic off- 
set in the A/, — Z relation. This is due to the fact that in our sim- 
ulation the probability to host a LGRB is higher in galaxies with 
higher SFRs. 

To further investigate these findings, we consider the FMR 
of the simulated catalogues dividing galaxies in subsamples with 
SFRs in different bins. The results are shown in the left panels of 
Fig lAll of Appendix A. We find that the FMR is well reproduced 
again by all subsamples of the HOSTl group, and by the subsample 
of galaxies with high SFR (log(SFR) > 0.9) of the H0ST2 group. 
By contrast, none of the HOST3 subsamples reproduce the data 
since the constraint on metallicity in the progenitors .Zprog implies 
also a severe cut in the metallicity of the host .Zhost . Therefore, we 
can conclude that collapsar model predicting a strong metallicity 
bias can not be easily reconciled with the observational evidence 
that LGRB hosts do follow the FMR. 



4 PROBABILITY DISTRIBUTIONS AND DARK LGRBS 

It has been noted that observed LGRBs tend to o ccur in galax- 
ies with high specific star f ormation rates (SSFRs) dSavaglio et al.l 
l2009l : lMannucci et alj20TTI) . The SSFR or its inverse, the so-called 
doubling time, are good marker s of the ability o f a galaxy to form 
stars in an efficient way. M annucci et al. N2OIII) have shown that 
the observed LGRB hosts have always doubling times shorter than 
the Hubble time at the redshift of the object. This holds true for 
additional 17 LGRB hosts at z < 1 of unknown metallicity for 
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Figure 2. Left: Mass-Metallicity relation for host simulated galaxies. The solid-black line is the Tremonti relation shifted to higher redshifts (Jabran et 
al.2010), the red-dashed line is the median for HOSTl, the yellow-long-dashed and green-dot-dashed line refer to H0ST2 and HOST3, respectively. Points 
are the observed host galaxies of our LGRB sample and the blue-dotted line is their best fit. Right: FMR for host simulated galaxies. We use the same line 
coding of the right plot. Black-solid line the FMR Mannucci et al. 2010a. 



which stellar mass M, and SFR have been measured 0. We note 
that none of them has log(SSFR) < -9.95 (with SSFR in units of 
yr~^) confirming that all observed LGRB hosts have high SSFRs 
and doubling times shorter than the Hubble time, at the redshift of 
the burst. The position of LGRB hosts in the SSFR-Z plane may 
thus provide new hints on the host galaxy population. 

In Fig. 3 we collect all observed data (blue dots) of the 18 
LGRB in the SSFR-Z plane (the 17 SSFRs for the additional hosts 
with lack of Z measurement are plotted as black point at the top 
of the figure). In the same plane, the color-coded area gives the 
normalized probability map to have a LGRB event in galaxies of 
the HOSTl sample (see Appendix for our analysis on HOSTl-2-3 
in separated bins of SFR or mass M,). The probability is computed 
normalizing the number of LGRB events of each galaxy for the 
total number of events. Most but not all observed LGRB hosts fall 
in the green zone where the expected probability for the HOSTl 
sample is the highest. The simulated galaxies cover a wide area in 
the SSFR-Z plane and extend also to low SSFRs. About 15% of 
all simulated LGRB s reside in galaxies with log (SSFR) < —9.93, 
i.e. with doubling times longer than the Hubble time at 2 = 0.5. 
Therefore the model does not exclude hosts with low SSFR and we 
may still lack of LGRBs due to the low statistics related to the bias 
against identifying galaxies with low SFRs, or to dust extinction 
(see below). 

Similarly, in the left panel of Fig. 4 we derive the color coded 
probability map of having a LGRB in a galaxy (with z ~ 0.5) of 
the HOSTl sample with given metallicity and /10.32 (the lack of 
simulated galaxies with /io.32 < 8.2 is due to the resolution limit 
of the N-body simulation used). Data-points are also shown for the 
observed LGRB host galaxies. As in Fig. 3, most of the simulated 
hosts cluster in a very well defined region of the FMR, and these 
loci are populated by the observed data. 



^ Data are from the GHOST database (www.grbhosts.org) 
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Figure 3. SSFR (in units of yr" ^ ) versus metallicity Z. Blue points are the 
data with associated e rror bars for the observed sample of 1 8 LGRB hosts 
iMannucci et alj201lL Black point at the top of the figure refer to 17 LGRBs 
for which the metallicity is not measured yet. The shaded area gives the map 
of the weighted probability to observe a GRB event for HOSTl sample. The 
lack of simulated galaxies with log SSFR > —8.5 is due to the resolution 
limit of the simulation used. 



We note the existence of a region (i.e., the green area in 
the interval 10 < /io.32 < 10.5 corresponding to 8.6 < 12 + 
log(0/H) < 9.1) with high values of probability of LGRB occur- 
rence that is underpopulated of observed LGRB hosts. The lack of 
observed hosts may be linked to the existence of dark LGRBs, i.e. 
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Figure 4. /io.32 — plane. In the left panel we show the color coded probability map of housing a LGRB for galaxies of HOSTl (with z ~ 0.5). Blue 
data points refers as in Fig. 3 to the observed LGRB hosts. In the right panel, we show the color coded probability map for a galaxy to have a given averaged 
extinction coefficient A^[mag] Simulated host galaxies are selected among the HOSTl catalogue at z ~ 0.5. The contours con'espond to the 50%(smaller 
area) and the 99% (bigger area) of the entire population. 



LGRBs with bright X-ray af terglow that are f aint or not detected 
in optical. It has been shown j Perlev et al.P 2009) that most of these 
objects can not be detected at optical wavelengths being strongly 
absorbed by dust in their host galaxies. Here, we argue that dark 
LGRBs may reside in hosts with high metal content (and therefore 
more dusty) that could populate the empty region in the ^0.32 — Z 
plane. Indeed the only point present in this region is relative to 
the host galaxy of ORB 020819, that is classified as a dark b urst 
jLevan et alj2003l : lKlose et alj2003l : [Kiipcu Yoldas et alj|2010l) . 

We briefly analyze the simulated galaxies (HOSTl) with high 
probability to host a LGRB in the empty region at high /io .32 , point- 
ing out that they have higher stellar mass (10 < log(M*/M0) < 
11), and they show higher extinctions. In the right panel of Fig. 4 
we show the color-coded map distribution of the total extinction 
coefficient A-o in the V band (Buser V filter). The contours contain 
the 50% (smaller area) and the 99% (bigger area) of the entire pop- 
ulation of simulated hosts showed in the left panel of Fig.|4] Again, 
the blue data points are the observed host galaxies. The principal re- 
sult of this figure can be summarized as follows: despite the higher 
probability for galaxies with high /^o.32 to house a LGRB, their 
higher extinction prevents us from observing the optical afterglows 
of their LGRBs, instead the galaxies with lower extinction, and thus 
lower metallicity, are favorable to be observed. The white and yel- 
low regions in the map refer to areas of low extinction where the 
bulk of the observed LGRBs are found. Similarly in the observation 
of core-collapse SNe, the missing in the Optical and IR observation 
is sup posed to be correlated with the extinction of dust in the host 
galaxylMaiolino et all ( l2002h : lMannucci et alj ( l2003h : [Cresci et all 
j2007h . 

However, we notice here that the present sample of observed 
LGRB hosts is not a complete or well controlled sample, and sev- 
eral selection effects could exist. In particular, we can not use the 
present sample to investigate the relative fraction of LGRB hosts 
detected as function of /io.32 Thanks to the unbiased simulated 
galaxy catalogues, we are able to predict the distribution on the 



FMR plane of host galaxies, waiting for future observations of dark 
GRB host galaxies. 



5 CONCLUSIONS 

In this paper, we studied the properties of LGRB host galaxies com- 
paring 18 observed hosts with a catalogue of simulated galaxies 
constructed combining high-resolution N-body simulations with an 
up-to-date semi-analytic model of galaxy formation. The simulated 
catalogue reproduces the offset on the observed M* — Z relation 
and the tighter Fundamental Metallicity Relation (FMR) of the field 
galaxy populations (in the redshift interval 0.07 < z ^ 0.3) with 
no need of introducing new free parameters. 

Aim of our investigation was at studying the relative impor- 
tance of metallicity and SFR in the characterization of a LGRB and 
of its host. To this purpose we studied the effects imposed by the 
presence of a metallicity threshold in the generation of a LGRB 
on the host galaxies. We extracted, from the main simulated cata- 
logue, three samples of LGRB hosts, housing different pockets of 
gas particles with different metallicities ^prog that could be poten- 
tial sites to LGRB formation. The first sample, denoted as HOSTl, 
comprises galaxies undergoing star formation but with no threshold 
on the metallicity (Zpio^). H0ST2 comprises instead star forming 
galaxies with lower metal content, i.e. a cut at Zprog = 0.3 Z©, 
and H0ST3 at Zp^og = 0.1 Zq. 

Our analysis shows that HOSTl is the sample providing a very 
good description of the currently available LGRB host dataset. In 
particular, the simulated sample is able to reproduce, at the same 
time, the systematic offset of LGRB hosts toward lower metallic- 
ities in the A/, — Z relation, and the tightness of the FMR in the 
/io.32 — Z plane. In the simulated sample, the probability of hous- 
ing a LGRB is higher in active star-forming galaxies. As a con- 
sequence, when weighted for this probability the simulated hosts 
tend to show lower metallicities at fixed stellar mass, in the M* — Z 
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plane. The HOSTl sample also reproduces the general trends in the 
SSFR-Z plane of observed LGRB hosts. By contrast, HOST2-3 
predict too lower metallicities in the A/* — Z relation with respect 
to what observed, and fail to reproduce the FMR of the available 
LGRB host sample. Also they exhibit an inverted behavior in the 
SSFR— Z plane with respect to the observed LGRB hosts. 

The close match of HOSTl with both the LGRB dataset and 
the FMR observed in the field, indicates that the request of a cut- 
off in the metallicity of the progenitor stars to LGRBs is not com- 
pelling. The low metallicity that is observed in the hosts of LGRBs 
appear to be a consequence of the higher probability of producing 
a LGRB in galaxies of low mass that have a higher star forma- 
tion rate. Since galaxies of a given stellar mass that are more star- 
forming are metal poorer, our analysis suggests the SFR to be the 
primary marker of LGRB formation. 

The observed sample of LGRB hosts is not complete and 
could be affected by strong bias effects. Indeed, the sample com- 
prises mainly LGRBs that displayed an optical afterglow. There ex- 
ists however, LGRBs with X-ray afterglows that do not show any 
optical counterpart, i.e. the so called dark LGRBs. Many of these 
are found in dust rich environments (Perley et al. 2009). It is unclear 
whether dust is spread over the galaxy or confined to the LGRB lo- 
cal environment. In the first case the host should have a high mean 
metallicity. These LGRBs should populate the loci of high metal- 
licity of the FMR that are devoid of objects, at present. The lack of 
such hosts is suggestive that the optical afterglow of such LGRBs in 
these galaxies suffer sever extinction, and so we are led to identify 
these missing events with the fraction of dark LGRB, according to 
our probability analysis study. 
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HOSTl follow the FMR. Galaxies with higher SFR and lower stel- 
lar masses have lower metallicities; also higher stellar masses cor- 
responds higher metallicities. 

In addition, as shown in Rig. lAll H0ST2 and H0ST3 sample 
provide a poor description of the distribution of observed LGRB 
hosts in the SSRR— Z plane. In these cases, the simulated galaxies 
show a flat or inverted behavior with respect to the data. Indeed, the 
metallicities of simulated hosts are found to increase with SSFR. 
This fact provide further evidence that models requiring a metallic- 
ity threshold for LGRB formation seems at odd with the properties 
of the observed LGRB host sample. 



APPENDIX A: DETAILS ON THE HOSTS' SAMPLES 

In this Appendix, we illustrate further findings obtained from the 
analysis of HOST 1-2-3 in selected SRR and A/, intervals. Rigure 
lAII show in the /io.32 — Z and SSFR— ^ plane the relations for our 
subsamples of host galaxies. The contribution of hosts with dif- 
ferent SRRs are shown with different colors. The same color cod- 
ing is adopted for the observed data. Rinally, the s hort-long dashed 
line is the best-fit to the observed data obtained by Mann ucci et al] 

1]). The FMR is well reproduced by all subsamples of galax- 
ies belonging to the HOSTl group, and they are also in agreement 
with observations (top panel of Fig. lAlt . For the H0ST2 sample 
(middle panel of Fig. lAlb simulated galaxies seem to provide a 
reasonable description of the data; both highly star forming objects 
(blue-dot-dashed lines and data) and less active galaxies (yellow- 
short-dashed and red-dotted lines) and data are well reproduced. 
However this model predicts a systematic shift of the LGRB host 
population with < log(SFR) < 0.9 (green-dashed lines) to- 
ward low metallicities with respect to the FMR. This shift is not 
observed in the present sample of LGRB host data (green points) 
that represent the bulk of the LGRB host population in the observed 
sample. This is reflected in the behavior of the total FMR of sim- 
ulated galaxies (solid line) that falls below the observed FMR of 
LGRB hosts. This effect is amplified in the H0ST3 sample (bot- 
tom panel of Fig. lAU Left Column), where only a few highly star 
forming galaxies with high metallicities are present. The shift in 
metallicity for middle and low star forming objects is more evident 
with the simulated curves always below the observed data points. 

Fig. lAII shows the distribution of simulated LGRB hosts at 
z ~ 0.5 in the SSFR— Z (right panels). The lines refer to the 
median weighted relations in different stellar mass bins, the data 
and lines of the same color refer to the same A/* range. Again, 
we find that HOSTl provides a good description of the data. In 
particular, HOSTl galaxies are found to have decreasing metallic- 
ities for increasing SSFRs and decreasing stellar masses similarly 
to the observed trends. This is a direct consequence of the fact that 
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Figure Al. Left column: /io.32 — Z plane for the host simulated galaxies (with z ~ 0.5), the colored lines are the median values in different bins of SFR. 
Data points refer as in Fig. 4 to the observed LGRB hosts and black-dashed line is their best fit. Finally black-solid line is the FMR Mannucci et al. 2010a. 
Right column: SSFR (in units of yr~^) versus metallicity Z. Points are the data as in the right panels. Triangles at the top of the figure refer to 17 LGRBs for 
which the metallicity is not measured yet. The lines are the median values for HOST 1-2-3 in different bin of Af» . 
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